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Abstract. Gamma-ray emission from AGN provides us with unprecedented insights 
into the physics of extragalactic jets. The emission from these jets fits naturally into 
the existing theoretical framework of relativistic jets as inferred from parsec scale and 
kiloparsec scale observations. Models of the 7-ray fluxes give us important knowledge 
J> ' of jet parameters, in particular magnetic fields and the size of the emitting region. 

. For example, in the case of Markarian 501, the jet plasma is subequipartiion and 

the energy flux is particle dominated. Discussion of the physics of jets on the sub- 
parsec scale leads naturally to questions of jet composition. The dynamics of jets 
are consistent with conventional electron-proton composition if the minimum Lorentz 
factors, 7 m i n > 100. Jets can be electron-positron if 7 min « 100. 
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O : INTRODUCTION 

CO 

One of the major discoveries of the EGRET experiment on the Compton Gamma 
Ray Observatory was the detection of significant fluxes of GeV 7-rays from radio 
loud active galactic nuclei (AGN). Since then ground-based TeV observations from 
the Whipple observatory and the High Energy Gamma-Ray Array (HEGRA) have 
ushered in a new era of ground-based 7-ray astronomy at the highest energies. 
Apart from one object - Centaurus A - all of the 7 rays detected so far from AGN 
have been from blazars. These are highly variable AGN resulting from jets oriented 
close to the line of sight. The aim of this talk is to give an overview of the physics of 
these objects, especially as they relate to high energy emission and the implications 
for processes involving relativistic jets. A specific TeV source, Markarian 501, is 
discussed in some detail. Much of this discussion leads naturally to the composition 
of extragalactic jets - electron-proton or electron-positron and some recent work on 
this question is also summarized here. 



^ To be published in Proceedings of the International Symposium on High Energy Gamma-Ray 
Astronomy, ed. F. Aharoninan and H. Voelk, American Institute of Physics. 



THE PARENT OBJECTS 



In the view of AGN known as Unified Schemes the parent, or unbeamed, coun- 
terparts of blazars are Fanaroff-Riley Class 1 and Class 2 radio galaxies (FR1 and 
FR2). FR1 radio galaxies contain initially relativistic jets which decelerate to ve- 
locities of 1, 000 — 10, 000 km s" 1 on the kiloparsec scale. Their properties on the 
kpc scale can be understood in terms of turbulent transonic flows. The deceleration 
from relativistic to sub-relativistic flow probably occurs on scales a few hundred 
parsecs. Consistent with this, FR1 radio sources are usually asymmetric close to 
the core but symmetric further out. This is interpreted as the diminishing effect of 
relativistic beaming. It is widely believed, and there is good evidence for this, that 
BL-Lac objects are FR1 radio galaxies viewed at small angles to the line of sight. 
However, sub BL-Lac objects may be beamed FR2s. 

The jets in FR2 radio galaxies appear to be relativistic along their entire lengths. 
Nevertheless, some deceleration is probably inevitable between the core and the 
hot-spots some tens of kiloparsecs distant from the core. The appearance of FR2 
jets is often dominated by knots - interpreted as internal shocks. They are almost 
always one-sided and Laing's and Garringtons's discovery [1,2] that the jetted side 
is almost always the least depolarized, provided a strong basis for the interpretation 
of their properties in terms of the beaming resulting from relativistic motion. Again 
in the Unified Scheme view of the physics of these sources, FR2 radio sources are 
viewed as the not highly beamed counterparts of quasars. 

COMPONENTS OF A THEORY OF 7 -RAY EMISSION 

FROM AGN 

The essential aspects of the mainstream models for production of 7-rays from 
AGN are the following: 

1. The unique association with radio-loud sources strongly indicates that the 
GeV-TeV emission is the result of inverse Compton scattering of soft photons 
by ultra-relativistic electrons in jets moving at relativistic bulk velocities. 

2. The inferred Lorentz factors of the relativistic electrons are as high as ~ 10 6 . 
Relativistic shocks (see Achterburg; Kirk and Mastichiadis, these proceedings) 
are the favored mechanism for acceleration of such high energy particles. 

3. A source of soft photons is required. These can be low energy photons from 
the jet itself (synchrotron self-Compton emission; SSC) or photons produced 
in the AGN environment (external inverse Compton; EIC) . (It is intriguing 
to reflect that here, "low energy" can mean X-rays which had at one stage 
represented the highest energy photons observed from AGN.) 

4. Relativistically moving plasma boosts the observed flux density by a factor of 
about 10 3 " 4 . 



Apart from models of the emission itself, there are other physical processes to 
take into account. These include: 

1. The high photon density required in SSC models implies a correspondingly 
high opacity related to pair production. This is especially important for TeV 
emission and is an important constraint on the size of the emitting region. 

2. What is the constitution of jets - electron-ion or electron-positron? 

3. What is the jet dynamics in the first 10 16 " 17 cm? Are jets initially Poynting 
flux or particle dominated? 

The Relative Importance of EIC and SSC 

Qualitatively, one may gain an idea of the relative importance of EIC and SSC 
emission from the energy density of the soft photons. Markarian 421 and 501 are 
blazars in which SSC models are favored. For those objects the radiation energy 
density is of order 10 -4 — 10~ 3 ergscm -3 . The radiation energy density at a distance 
10 17 ri7 cm from an accretion disk emitting a luminosity of 10 44 L44 ergs s" 1 is 

Mdisc ~ 0.03^- ergs cm 
r( 7 

This expression does not describe exactly when EIC or SCC emission may be 
important. For example, [3] showed that when scattered diffuse radiation from 
an accretion disc is taken into account, the radiation density is enhanced by a 
factor ~ T 2 . Nevertheless, there is little evidence in FRls and BL-Lac objects for 
significant accretion disk luminosity. Hence, it is not surprising that for BL-Lacs, 
SSC models are often favored. Moreover, Dr. Krawczynski has also presented 
evidence at this meeting that the TeV variability of Markarian 501 is quadratic in 
the X-ray flux, supporting an SSC model for that object. 

Synchrotron Self-Compton Radiation 

Historically, synchrotron self-Compton radiation was considered in the context of 
the production of X-ray radiation from radio-emitting plasma, the inverse Compton 
catastrophe and the limit of 10 12 K on the brightness temperature of a synchrotron 
emitting plasma [4]. SSC radiation at these energies is generally produced in the 
Thomson limit wherein the energy of the scattered photon ~ 7 2 e where 7 is the 
electron Lorentz factor and e is the energy of the soft photon. Thus, 7 ~ 10 4 
electrons produce 1 keV photons from GHz photons. We know now that higher 
energy photons, from infrared through to X-ray are produced from synchrotron 
emission at the sub-parsec - parsec bases of jets and the corresponding SSC photons 
are in the GeV to TeV range. In some case the emission can extend beyond the 
Thomson limit to the Klein-Nishina regime. For the latter case the energy of 



the scattered photons is of order ^m e c 2 . Allowing for a Doppler factor of 10, 
say, observed emission at 10 TeV implies a Lorentz factor of about 2 x 10 6 . It 
is generally envisaged that electrons with such high Lorentz factors are produced 
locally at shocks. The shocks in blazar jets are probably weak to moderate since 
it is not necessary to convert a large amount of the jet energy flux into emitting 
particles to account for the observed fluxes. 



Boosting of Emission from Relativist ically Moving Plasma 

The Doppler factor, S of an object moving with velocity v — c(5 at an angle 6 to 
the line of sight, is: 

s = —^ .. 

r(l -/3cos0) 

where the bulk Lorentz factor is Y = (1 — f3 2 )~ 1 ^ 2 . Let the observed flux density 
from a region of relativistically moving plasma be F„, the luminosity distance be 
Dl , the emissivity at the emitted frequency, v' be j v i and the element of volume 
in the rest frame be dV, then 




This is one of the most fundamentally important equations of blazar physics and 
clearly shows that for Doppler factors of 10 — 20 that are inferred in many blazars, 
there is a factor of 10 3-4 enhancement in flux density over that for a Doppler factor 
~ 1 source. This, of course, is the lynchpin of Unified Schemes which have been so 
important in linking the various members of the AGN zoo. However, it is attractive 
to think of this relationship in a different way and that is as a method of using 
the relativistic boosting of blazars to analyze the normally faint regions close to 
the black hole in radio-loud active galaxies. This is analogous to the way in which 
optical astronomers utilize the amplification resulting from gravitational lensing to 
observe galaxies at high redshift. It is interesting to comtemplate that both effects 
are consequences of the special and general theories of relativity. 



THE DYNAMICS OF RELATIVISTIC PLASMA 

Superluminal motion 

Observations at high resolution by astronomers on physical scales ranging from 
parsecs to kiloparsecs provide information that is relevant to processes occurring 
on the scales accessed by 7-ray astronomy. First, of course, observations of super- 
luminal motion have been recorded by VLBI radio astronomers for many quasars 
and BL-Lac objects. Moreoever, Biretta and colleagues have observed superluminal 



motion close to the nucleus of M87 using the HST [5]. The apparent superluminal 
velocity is of course linked to the actual velocity by the well known expression: 



(3 sin6 
' app ~ l-^costf 

with a maximum value of Tf3. Thus observations of superluminal motion provide 
a lower limit on the Lorentz factor. 



Pattern speed and jet speed 

The observation of the apparent transverse speed of a moving knot in a jet does 
not immediately translate to a constraint on the velocity of the jet. The moving 
knots are probably shocks [6] whose pattern speed is likely to be different from that 
of the jet. For example, consider a weak jet shock which, in the frame of the jet 
plasma, is moving towards the core. It is therefore swept out from the core by the 
jet but at a lower velocity than the jet. For a weak shock, take /3j e t, s h ~ l/\/3 and 

Fjet,sh ~ ^3/2 to be the velocity and Lorentz factor of the jet plasma relative to 
the shock, (3 sh xs 1 and r sh the laboratory frame shock velocity and Lorentz factor, 
then the transformation between the shock and laboratory frames tells us that the 
lab. frame Lorentz factor of the jet is 

Tjet = r s hFj etiSh [1 + /5 s h/9jet,sh] ~ 1-9 r sh 

The factor of 1.9 increases for stronger shocks and underlines the fact that even 
weak shocks have a substantial velocity in a relativistic plasma. 

For a jet observed almost pole-on (6 ss 0), 5 2Fj Ct 3.8r sh for a weak shock. 
Hence a Doppler factor of 20, as inferred for MKN 501 (see below) is not unexpected, 
given that superluminal velocities of the order of 5 c are common. 



LESSONS FROM KPC SCALE JETS 

The famous radio galaxy M87 is close enough that high resolution observations 
over a period of 11 years [7] have successfully detected motions of order 0.5c. More 
recent HST observations have revealed superluminal motions close to the core. A 
velocity of 0.5c is not fast enough to be consistent with the required sidedness 
and this motivated the interpretation of the bright knots in the M87 jet as reverse 
shocks moving more slowly than the underlying jet plasma [8] . The jet also oscillates 
significantly beyond knot A. Begelman and Bicknell interpreted this as a body mode 
Kelvin-Helmholtz instability. The pattern speed of the knot is the pattern speed of 
the instability and a plot of maximally growing wavelength against pattern speed 
is shown for a T = 5 jet in figure (1). The range of values indicated by the data is 
shown. The interesting point to note here is that the best fit between theory and 
data is only attainable when 



Rest mass energy density pc 2 

X = = — ~ o — 1 ( 1 ) 

Enthalpy 4p 

This is only possible if the jet consists of electron-positron pairs or if it consists of 
electron-proton plasma with a relatively high lower cutoff in the electron Lorentz 
factor. Specifically, for a relativistic electron distribution described by ^(7) = 
K e ^~ a for 7 min < 7 < 7 m ax, the value of x for an electron-proton plasma is 

3 m n (a — 2\ , 
Xep = j— ( — — r 7min « 230^ for a = 2.2. (2) 

Thus for M87 and electron-proton jet requires 7 m i n ~ 200. On the other hand a 
simple alternative is that these jets are composed of electron-positron pairs. We 
address this issue further below. 



LESSONS FROM MARKARIAN 501 



This section is the result of research carried out in collaboration with Stefan 
Wagner (Bicknell & Wagner, in preparation). 



X-Ray and 7-ray light curves of Markarian 501 

Markarian 501 is one of the two confirmed TeV blazars, the other of course, 
being Markarian 421 [9]. Both of these objects are good examples of high energy 




FIGURE 1. Plot of the wavelength of the maximally growing Kelvin- Hclmholtz mode against 
pattern speed of the instability for jets with varying values of the parameter \ = pc 2 /Ap and values 
of r = 3,5. The point and dashed line indicates the range of values implied by the observations 
of the M87 jet. 



SSC emission from blazar jets. Radio images [10] show a typical one-sided core- 
jet structure with large apparent jet deflection consistent with an almost pole-on 
jet. Two epoch observations have indicated a superluminal velocity of 6.7c for one 
component. However, this requires confirmation by observations at more epochs. 

We have been undertaking analysis of multi-epoch observations of MKN 501 over 
a period where it was flaring. The combined RXTE and HEGRA observations show 
an intriguing delay between the 1 keV X-ray minimum and the TeV minimum. 
(The latter occurs approximately one day later.) The reason for this behaviour 
becomes apparent when one examines the evolution of the X-ray spectra above 
10 keV. These show substantial hardening at the same time as the 1 keV emission 
is diminishing. Clearly, these observations caught MKN 501 at a time when one 
flare was subsiding and another was starting to develop. Unfortunately, as the 
second flare was increasing, the moon was rising and the flare could not be tracked. 
A detailed analysis of this data will soon be published. However, some of the 
implications for the parameters of MKN 501 are already apparent. 

Constraints on parameters of emitting region 

Let us consider some of the constraints on the parameters of the emitting region 
which are common to many analyses of blazars [11,12,9] in the context of a spherical 
geometry for the emitting region. 



There are two features of this equation. First, even for relatively large Doppler 
factors, the size of the emitting region is quite small, an order of magnitude smaller, 
for example, than the size of knots regularly observed with VLBI. Secondly, it is 
large Doppler factors that save those emitting regions which vary on a time scale 
of less than day, from becoming too compact and pair opaque. 



As Mattox [13] pointed out, the production of pairs by photon interactions (pair 
opacity) is important for compact 7-ray emitting regions. As Renault has also 
shown at this meeting (see also [11]), pair opacity of 7-rays on the IR background 
is also important. 



Variability 



A variability time scale of At implies that the radius, 





Pair opacity 



The threshold energy of soft photons for pair production on 7-rays with an ob- 
served energy, e 7 , is: 



Csoft 



2 4 

m e c 



0.265 



TeV 



The photon-photon cross-section peaks at the threshold energy so that a substantial 
photon density at the threshold energy can be important for 7-ray absorption. 
Using expressions calculated by [14] for a power-law electron distribution, the pair 
opacity for a spherical emitting region is: 



t « 2.6 x l(r 



F„(keV) 
/iJy 



i?r6^" (2Q+3) z 2 



TeV 



0.5 



It is implicit in this equation that the observed X-ray flux density is used to ex- 
trapolate the photon density to the infrared. It is the high power of the Doppler 
factor in this expression which rescues TeV blazars from potential obscurity. Note 
also the dependence on the radius of the emitting region, implying a lower limit on 
R, for the detection of 7-rays. Emission at TeV energies, of course, also provides 
much more severe constraints on size than emission at GeV energies. 



Spectral Breaks 

Both MKN 421 and 501 exhibit spectral breaks in the X-ray. In MKN 501 the 
break in spectral index is Aa 0.3 at about 4 keV. Such a break is probably 
the result of post-shock cooling and this aspect of blazar physics is not strictly 
consistent with the assumption of a completely homogeneous emitting region. The 
region following the shock must be stratified. For an homogeneous region of finite 
extent, it can be shown that the superposition of cooling spectra imply Aa = 0.5. 
In an expanding region following a shock, the contribution of the low energies 
to the integrated spectrum is diminished by adiabatic cooling and one can see, 
qualitatively, that the magnitude of the spectral break is decreased. Models of this 
effect, could further constrain blazar parameters. 

The location in energy of the spectral break can be estimated using the criterion: 

Travel time across emitting shocked region = Cooling time 

The travel time is estimated in the frame of the shock and one takes into account 
the time dilation of plasma moving in this frame since the synchrotron cooling time, 
for example, is estimated from the magnetic field in the plasma rest frame. If we 
take account of synchrotron cooling only, then, the radius of the spherical emitting 
region is given by: 



where, (f3 f , T') are the velocity and Lorentz factor of the plasma in the shock frame. 
The quantity, T'f3' m 0.71 for a weak shock and 0.35 for a strong shock. The 
numerical value is for a weak shock. Usually, the inverse Compton cooling time is 
comparable to the cooling time so that the cooling length is approximately a factor 
of two shorter. 



An SSC fit to X-ray and TeV epochs of Markarian 501 at which the emission 
is likely to be dominated by the rising flare is shown in figure 2. In this fit the 
TeV data have not been corrected for intergalactic pair opacity. This increases 
the magnetic field and decreases the size of the emitting region but not by a large 
factor. 

As you can see from the fits to the MKN data two components are required to 
model the X-ray part of the spectrum consistent with our earlier remarks concern- 
ing the emergence of a hard component. The cutoff Lorentz factor, 72, for the 
hard component is determined by the peak of TeV spectrum and the fact that the 
emission is in the Klein-Nishina limit. Hence the peak of the TeV spectrum is at 
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FIGURE 2. A fit to corresponding epochs of the X-ray and TeV data for Markarian 501 for 
a Doppler factor of 10. Two components are required to account for the hard and soft X-ray 
components. 



an energy of 57 2 m e c 2 . In this fit 72 = 3 x 10 6 ; the magnetic field is 0.06 G and the 
radius of the emitting region is 1.3 x 10 16 cm. 

The upper limit on the radius inferred from a variability time scale of one day 
(see equation (3)) is exactly the same as the radius inferred from the model. The 
cooling radius estimated from equation (4) above ~ 7 x 10 15 cm. Better consistency 
with the constraints imposed by cooling can be obtained with a Doppler factor 
~ 20 and the upper limit on R implied by variability is more comfortably satisfied. 
However, these parameters are probably in the right ballpark and are similar to 
those obtained by [12] for an earlier flare. 



PHYSICAL INFERENCES FROM BLAZAR MODELS 

The Milano group (Celotti, Fossati, Ghisellini, Maraschi and collaborators) have 
also obtained fits to numerous blazars in a similar fashion. The parameters obtained 
can give us valuable insight into the physics of jets on the sub-parsec scale and again 
we use Markarian 501 as an example. 



Relative energy of particles and magnetic fields 

The magnetic energy density of the emitting region in Markarian 501 is ap- 
proximately 1.6 x 10~ 4 ergs cm" 3 and the particle energy density is approximately 
3.7 x 10~ 2 ergs cm -3 . Thus the jet is subequipartition by a large factor. This has 
ramifications to models of jets [15] in which magnetic processes drive the jets and 
produce a flow in which there are equal amounts of Poynting flux and particle flux. 
Either issues of particle loading have to be taken into consideration or different 
mechanisms for jet launching have to be considered. 



Energy flux and electron-positron jets 
The energy flux of a jet consisting of relativistic particles plus cold matter is: 



1 E 
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where A = irR 2 is the cross-sectional area of the jet and as above x — pc 2 /Ap is the 
ratio of cold matter rest energy density to enthalpy [16]. For Markarian 501, 



7.5 x 10 43 
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For x ~ this is typical of (actually slightly greater than) the energy that we would 
associate with a BL-Lac object or (assuming unification) an FR1 radio galaxy. 
Independently then, we again arrive at the option of either an electron/positron jet 
or an electron/proton jet with 7 min > 100. 



The Gammasphere 



The notion of a gammasphere was introduced by Levinson and Blandford [17] in 
a theory for the EGRET sources. For the parameters relevant to Markarian 501, 
the pair opacity is 



The emitting region in Markarian is only just transparent for a Doppler factor of 
10 and quite transparent for a Doppler factor of 20 (r ~ 0.05) again indicating a 
Doppler factor of this magnitude. Thus the TeV 7-ray emitting region in Markarian 
501 is not too far above the distance from the black hole where it would become 
opaque. This raises the question, for Markarian 501, and for other blazars, as to 
whether there are X-ray flares which are not accompanied by 7-rays. 



At two places in the above description of the physics of relativistic jets the ques- 
tion of jet composition has arisen. It is therefore worth considering this question in 
some detail. The main issue to be addressed here is whether an electron-positron 
jet can be powerful enough to supply the requisite energy to the outer parts of a 
radio galaxy, quasar or BL-Lac object. An important physical constraint is that 
annihilation limits the jet density and this may be inconsistent with what is re- 
quired by the energy budget [18,19]. Celotti and Fabian cleverly worked out a 
way to independently estimate the energy flux of a jet by exploiting a previously 
discovered correlation between jet energy flux and narrow line region luminosity 
[20]. Groves and Bicknell (in preparation) have revisited this question essentially 
by adding extra terms to the energy flux omitted by Celotti and Fabian. 

In our sample of quasars and BL-Lac objects (mainly the sample of Celotti and 
Fabian with a few additions) the inverse Compton emission from these sources is 
assumed to be X-ray emission and the soft photons are from the radio part of 
the spectrum. (This led to some interesting discrepancies when some well known 
X-ray-7-ray blazars crept into the sample!) 

The parameters of these objects were estimated following the approaches in 
[18,19], that is: 

• The ratio of X-ray to radio flux density is used to estimate the electorn density 
(specifically the parameter K e in N(^y) = K e ^y~ a ). 

• The synchrotron self-absorption turnover, the synchrotron flux and the VLBI 
angular size are used to estimate the magnetic field and the Doppler factor. 

• The minimum Lorentz factor was estimated from the Doppler factor. 




ELECTRON-POSITRON JETS 
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FIGURE 3. Jet energy flux plotted against narrow line region luminosity for the Rawlings and 
Saunders radio galaxies (filled circles) and for the quasars and BL-Lac objects in the expanded 
Celotti and Fabian sample. 



For the present purposes we express the energy flux as: 

F E = c (3A ]ct [(r 2 -r) pc 2 + 4r 2 p 

The first term represents the kinetic energy flux associated with the rest-mass den- 
sity (of either cold matter or the relativistic particles); the second term represents 
the enthalpy flux associated with the relativistic particles. This term can be the 
most important term in the energy flux. Take n e to be the density of both elec- 
trons and positrons then pc 2 = n e m e c 2 and, as before, take Nfr) = K e ^y~ a for 
7min < 7 < 7max to describe the distribution of Lorentz factor of the relativistic 
electrons and positrons. Then, the relative importance of enthalpy and kinetic 
fluxes is given by 

4p _ 4 a - 1 

n 77 7r7min ~ o7min 

nm e c z 6 a — 2 

for a = 2.2. Hence, for large values of 7 m i n the enthalpy term dominates and we 
are able to estimate a large energy flux for a modest density. Interestingly, the 
viability of an electron-positron jet model depends on a large value of 7 min as does 
the viability of an electron-proton model. 
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FIGURE 4. Electron-positron density plotted against the survival density of electrons and 
positrons. Filled circles correspond to 7 m ; n = 10 and open circles correspond to 7 m j n = 100. 

Figure 3 shows the plot of jet energy flux against narrow-line luminosity for 
7min = 100. The energy fluxes in our sample extend this correlation quite nicely 
and this results from the importance of the enthalpy flux. 

The crucial plot is that of pair density. Celotti and Fabian used an estimate of 
the number density of surviving pairs at the VLBI scale. Whilst this estimate is 
model-dependent it is probably quite good since the main effect is the competition 
between annihilation and expansion. We have used the same model and in figure 4 
we plot the number density of electrons and positrons against the survival number 
density, for 7 min = 10 and 100. For 7 min = 100 the number density is comfortably 
below the survival density. These two plots therefore show that an electron-positron 
jet is therefore viable for this value of 7 m i n . There is a distinction here from the 
electron-proton case which requires r y m in > 100. 

DISCUSSION 

The gamma-ray emission from active galactic nuclei fits comfortably into the 
picture of relativistic jets that has been developed from both a theoretical and ob- 
servational perspective, mainly inspired by radio observations, over the last twenty 
five years or so. An appealing aspect of blazars is that, through the effects of 



relativistic beaming, they give us a way of investigating the physics of relativistic 
jets close to the environs of black holes and lead us to confront some of the major 
issues facing black hole physics in active galactic nuclei. These are: How are jets 
produced; what are the physics of jet propagation in the first 10 16 cm and why 
do some AGN produce jets and why do some not? The developing field of 7-ray 
astronomy, one of the last frontiers, is well placed to provide important information 
on these issues. One looks forward to the next generation of 7-ray telescopes and 
to attempting to resolve many of the problems with which they will surely provide 
us. 

GVB is grateful to the organisers of this meeting for assistance in funding his 
visit to Heidelberg. 
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